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£T) I ABSTRACT 

■ Using deep, high-spatial resolution imaging from the HST ACS Coma Cluster Trea- 

sury Survey, we determine colour profiles of early-type galaxies in the Coma cluster. 
From 176 galaxies brighter than M F814W (ab) = — 15 mag that are either spectroscop- 
ically confirmed members of Coma or identified by eye as likely members from their 
low surface brightness, data are provided for 142 early-type galaxies. Typically, colour 
. profiles arc linear against log(i?), sometimes with a nuclear region of distinct, often 

J-j ' bluer colour associated with nuclear clusters. Colour gradients are determined for the 

regions outside the nuclear components. We find that almost all colour gradients are 
negative, both for elliptical and lenticular galaxies. Most likely, earlier studies that 
report positive colour gradients in dwarf galaxies are affected by the bluer colours 
of the nuclear clusters, underlining that high resolution data are essential to disen- 
tangle the colour properties of the different morphological components in galaxies. 
Colour gradients of dwarf galaxies form a continuous sequence with those of elliptical 
galaxies, becoming shallower toward fainter magnitudes. Interpreting the colours as 
metallicity tracers, our data suggest that dwarfs as well as giant early-type galaxies in 
the Coma cluster are less metal rich in their outer parts. We do not find evidence for 
environmental influence on the gradients, although we note that most of our galaxies 
are found in the central regions of the cluster. For a subset of galaxies with known 
morphological types, SO galaxies have less steep gradients than elliptical galaxies. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: dwarf - galaxies: clus- 
ters: individual: coma - galaxies: structure 
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Metallicity gradients provide a means of studying galaxy 
formation. Classical monolithic collapse scenarios (|Larsonl 
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ll974l ; ICarlb"er3ll984l ; lArimoto fc Yoshiill 1987ft predict strong 
metallicity gradients. In this scenario primordial clouds of 
gas sink to the centre of an overdensity where a rapid burst 
of star formation occurs. Infalling gas mixes with enriched 
material freed from stars by stellar evolutionary processes 
and forms a more metal rich population. Because the gas 
clearing time (and hence the number of generations which 
enrich the interstellar medium) is dependent on the depth 
of the potential well, the metallicity gradient is dependent 
on the mass of the galaxy. 

Mergers, dominant in hierarchical galaxy formati on sce- 
narios, will dilute existing p opulation gradients (e.g. IWhitel 
1 198(1 Id i Matteo et al J 120091 ). although resid ual central star 
form ation can steepen gradients again (e.g. iHopkins et al.l 
2009, among others) . Therefore, the study of metallicity 
gradients can be used to distinguish between competing sce- 
narios of galaxy formation and can eventually lead to a more 
detailed understanding of those scenarios. Many attempts 
have been made to model metallicity grad i ents in more detail 
ilCarlberg||l984l; lArimoto fc Yosmt Il987l; jChiosi fc Carrarol 
120021 ; iKobavashil |2004 iPipino et all |200S| , l2010ft . of which 
the models bv lKawata fc Gibsonl (120031) are so far most con- 
sistent with observations (ISpolaor et alj|2009ft . 

Much work has been devoted to determine metal- 
licity gradients in lar ge elliptical gal a xies, either b y use 
of colour gradient s llSandagel Il972l; iFranx et al.l 19891; 
Peletier et~atl Il990ll?l: ISaglia et all I2OO0I; lLa Barbera et al.1 



20041 . 120051: lLa Barbera fc de Carvalhol 12009ft or by spec- 
troscopy (iMehlert et aIll2003l:ISanchez-Blazauez et al-lboOfl 
iQgando et al.ll2008l ; ISpolaor et alj|2010l ). There is consensus 
among all authors that elliptical galaxies do have metallic- 
ity gradients, which are not as strong as predicted by clas- 
sic monolithic collapse scenarios. However, for dwarf g alax- 
ies, the situation is not as clear. IVader et al.1 (|l988ft find 
that colour g r adien ts become positive for dwarf galaxies. 
Ivan Zee et al.l (|2004l ) mitigate this statement by saying that 
non-nucleated dwarf galaxies have positive gradients, but 
nucleated dwarf ga laxies have more or less flat gradients. 
iGorgas et al.l |l997ft find negative metallicity gradients for 
a sample of high-luminosi t y dw arf galaxies. The metallicity 
gradients of lSpolaor et al.l (|2009l ') show a trend to flat or pos- 
itive gradients at low luminosities, whereas the metallicity 
gradients from iKoleva et all (|2009) stay again negative in 
this regime. 

In addition to the different formation scenarios, differ- 
ent evolutionary paths may have played a role in shaping the 
dwarf galaxies in different parts of the cluster. It is obvious 
that at some point dEs have started as entities containing 
gas, however, as we observe th em now they ar e gas-poor. The 
morphology-density relation ( Dressie r 1980) suggests that 
cluster-specific processes are important for removing the gas. 
Several different processes have been proposed to operate in 
clusters (and less dense environments), but it is not clear 
which one is the most imp ortant in the end. Ra m pressure 
stripping by the hot ICM (|Gunn fc Gottl 1 1972ft will result 
in a differ ent stellar po pulation than f or example harass- 
ment (| Moore et al.lll99q ) and starvation (|Tullv fc Trent haml 
2008). By looking at gradients in different environments it 
may be possible to identify the most promine nt evolutionary 
paths (if such a thing exists, see for example l|van Zee et al.l 
12004ft .) It is obvious that this is extremely difficult, because 
all gas removing processes act more efficiently on gas in the 



outskirts of the galaxies and are more inefficient in the re- 
moval of gas in the centre. If the left-over gas in the centre 
forms stars, this will steepen existing gradients a bit, but it is 
not clear by how much (except in the case of tidal stripping, 
which removes also part of the stellar content of a galaxy) . 
However, it is clear that these processes depend on galaxy 
density or position in the cluster. The Coma cluster is there- 
fore an ideal laboratory to study the evolution of galaxies, 
as the density is to a first approximation only a function of 
clustercentric distance. 

Coma, at a distance of 100 Mpc the nearest Abell class 
3 cluster, is an excellent place to study dwarf galaxies for 
various reasons. Since Coma is located at high G alactic lat- 
itude, the extinction is low (Sc hlegel et al. 1998). Further- 
more, Coma, serving as a benchmark for high redshift stud- 
ies, has been looked at intensively in the past - both with 
ground based telescopes and with space-based telescopes 
(GALEX, Spitzer) - and in the present. With the advent of 
the Coma Treasury survey many new observing programs 
have commenced|j 

In this paper we study colour gradients for dwarf and gi- 
ant galaxies in the Coma cluster as derived from HST data. 
Colours probe stellar populations: colours of metal rich stars 
become redder, but also the age of the population and the 
slope of the initial mass function (IMF) change the colour. 
The superior resolution of HST allows us to separate differ- 
ent morphological components of galaxies, leading to a clean 
interpretation of our observations. 

Throughout this paper we assume a distan ce modulus 
to Coma of m — M = 35 |Carter et al.l 120081 . henceforth 
paper I), which corresponds to a distance of lOOMpc. 



2 DATA AND OBSERVATIONS 

2.1 HST/ ACS Coma Cluster Survey data 

We have used the data from the HST/ ACS Coma Cluster 
Survey, a deep 2-passband imaging survey of the Coma Clus- 
ter. A full description of the observations and data reduction 
can be found in Paper I. Here we will give only a short sum- 
mary of the survey. The ACS Wide Field Camera was used 
to image the Coma Cluster in the F814W and F475W pass- 
bands (roughly coinciding with the I and g filters) . The sur- 
vey was designed to contain different galaxy density regimes, 
by observing both the core of the cluster and the outskirts, 
and in particular the infalling group of NGC 4839. Due to 
the breakdown of the ACS Camera the survey is only partly 
completed; with 19/25 observed fields the centre of the clus- 
ter is well-covered. The other six fields sample the region 
around the infalling group NGC4839. Figure [T] shows the 
positions of the frames on the sky. The fields that were com- 
pletely observed (4/25 fields are only partly observed) had 
exposure times of 4 x 350s in F814W and 4 x 640s in F475W. 
After data reduction, which included PyRAF/STSDAS and 
MultiDrizzle (see paper I for details), each ACS frame con- 
sists of 4300x4225 pixels, each having a size of 0.05" xO.05" . 
A source catalogue was created by running SExtractor on 



1 An overview of observing programs can be found at 
http : //astronomy. swin.edu. au/ coma 
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Figure 1. Left: Positions of the Coma ACS fields, overlayed on an XDSS image. The length of the arrow of the compass is 10'. Right: 
position of sample galaxies. 



the reduced frames. A fu ll description of the source cata- 
logue is given in paper II (|Hammer et al-lfeuTOt ) 



2.2 Sample selection 

Except for the galaxies in the most central fields, the major- 
ity of observed galaxies are background galaxies that do not 
belong to the cluster. For our analysis we are only interested 
in galaxies inside the Coma Cluster. Arguably the most reli- 
able way of determining cluster membership is by determin- 
ing the redshift of each galaxy. For a large sample of galax- 
ies this has been done using the Hectospec sp ectograph on 
MMT , and t he LRIS spectrograph o n Keck fsee lSmith et al.l 
l|200Sl . 120091 ) : IChiboucas et all (|2010l ').') The sample for which 
we have accurate redshift information is approximately com- 
plete for galaxies with /i e (F814W) < 22.5. 

The majority of the fainter Coma dwarf galaxies stand 
out against background galaxies as low surface brightness 
galaxies. Due to the superior resolution and sensitivity of 
the ACS camera it is therefore possible to estimate clus- 
ter membership solely on the basis of morphology for faint 
galaxies of which membership cannot be confirmed spectro- 
scopically. We use the same sample as used for the determi- 
nation of the luminosity function for Coma. For details on 
the selection method and the selected galaxies, see Trentham 
et al. (2010, in preparation). The sources selected in Tren- 
tham et al. all have a membership rating, going from to 3. 
Class sources are spectroscopically confirmed members of 
Coma. Class 1 contains "almost certain" members, which are 
sources that are only rarely found in the field. Class 2 con- 
tains "likely" members, and class 3 contains possible Coma 
members. Spectroscopic follow-up has shown that rating 
class 2 contains approximately 10% non-members, whereas 
rating class 3 contains ~50% non-members. To minimise 
background contamination, we only use members from class 
(confirmed) to class 2 (likely member), but leave out class 
3. Our sample contains approximately 176 galaxies, 10% of 



which are of class 2. This means that the expected number 
of non-members in our sample is 2. 

Not all galaxies in this sample can be used for determin- 
ing colour gradients. First of all, we are only interested in 
early-type galaxies, and we therefore exclude galaxies with 
obvious spiral features. For sources fainter than F814W(AB) 
= 20 the dwarfs become generally too faint to obtain any 
reliable results at all, so we also have excluded these sources. 
In addition to this, the sample contained a number of ex- 
tremely compact galaxies (UCDs) which were picked up by 
Hectospec and LRIS. Since their size is hardly larger than 
the size of the PSF, we cannot obtain any useful colour gra- 
dients from them. We have excluded galaxies where the pres- 
ence of dust in the outer parts of the galaxies prevented us 
from probing stellar population differences. In addition to 
this, there are a number of galaxies for which it is impossi- 
ble to obtain reliable results, because they reside too close 
to other galaxies, and it is impossible to get a reliable es- 
timate of the sky, or because the shape of the galaxy is 
such (highly flattened, or irregular) that the fitting proce- 
dure does not give very accurate results. Table IA1I in the 
Appendix contains the list of 142 fitted sources and their 
membership rating. Table IA2I contains a list of 36 galaxies 
that have been excluded for various reasons. Out of these 
142 sources, 1 source has membership class 1, and 16 have 
membership class 2, meaning that we have cluster member- 
ship confirmation for 125 out of 142 sources (88%). 



3 ANALYSIS 

3.1 Galaxy isophotal analysis 

We use Galphot (|Franx et aHll989l ; |j0rgensen et ai1ll992l ) to 
measure the surface brightness profile of each galaxy. This 
program fits ellipses to a light distribution, allowing a vari- 
able centre, ellipticity and position angle as function of ra- 
dius. 
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For our colour gradients, all the fitted ellipses need to 
have a common centre. The centre as determined by Sex- 
tractor is in many cases not accurate enough, because its 
luminosity weighted determination will not coincide with 
the real centre, for example, if a galaxy is slightly asym- 
metric or if the sky background has some gradient. In this 
case one can either choose to centre the ellipses, after an 
initial Galphot run, on the outer isophotes or on the inner 
isophotes, or some combination of both. Motivated by the 
presence of nuclei in most dwarf ellipticals, we decide to 
fix the centre on the inner ellipse. We note however that a 
slightly different centre has only a negligible effect on the 
gradient. 

The shape of the central light profile of a galaxy is usu- 
ally quite sensitive to the size of the point spread function 
(PSF). To obtain useful colour profiles one has to compare 
fluxes at the same place in the galaxies. We convolved each 
image stamp in each ba nd with the psf of the other bands. 
PSFs were generated fcf. lHovos et alj (|201Ch . henceforth Pa- 
per III) using Drizzly TirrQ. Convolution was done in fre- 
quency space using a Fast Fourier Transform. 

In summary, we set up our Galphot runs in the fol- 
lowing way. In the first run, we give Galphot as input the 
centre of the source as determined by SExtractor, and al- 
low the centres of ellipses to change. This Galphot run is 
done on a convolved F814W image stamp. In the second 
run, we use the results from the first run to fix the centre, 
ellipticy and position angle, such that they are the same for 
all ellipses. After this, we determine the azimutally averaged 
surface brightness of the same galaxy in the F475W band, 
using the same set of ellipses. A colour gradient is then deter- 
mined by fitting a straight line in colour-log(r) space, where 
r is the circularized distance of the ellipse (i.e. the geometric 
mean of the major and minor axis.) This scale-free way of 
fitting avoids the use of a predetermined scale in a galaxy, 
such as the effective radius, and is common in the literature. 
Particularly for galaxies with strong gradients using a loga- 
ritmic gradient gives a better fit than a linear gradient. The 
range over which we fit extends from the innermost point to 
the cutoff in the profile (see section EO)) . and the points are 
weighted in a standard fashion by the inverse of their stan- 
dard deviation (a combination of the error in the sky and 
the root-mean-square deviation inside each elliptical annu- 
lus). The errors on the gradients quoted in this paper are 
the formal errors from the fit. 



3.2 Sky determination and uncertainty 

Arguably the hardest part of the analysis is the determi- 
nation of the sky background. In our observations, the sky 
background is made up of zodiacal light and earth shine, but 
on top of that, diffuse intra-cluster light, tidal debris, scat- 
tered light from bright objects and envelopes of big galaxies 
affect the measured background, so that determining the sky 
is an extremely difficult problem. 

SEXTRACTOR provides a measure of the sky and also the 
standard deviation of the sky. Nevertheless, it is well known 
that the apertures used by SEXTRACTOR are in general too 
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small for large galaxies (see e.g. Paper III), leading to an 
overestimation of the sky. 

We determine the sky value in the following way. The 
host galaxy and all other sources in the cutout are masked. 
We randomly pick 9x9 pixels sized boxes in non-masked 
areas, for which we determine the median. The median of 
the boxes is then used as the sky value. We find that the 
uncertainty in the sky is typically 2% in the F814W (and 
smaller in F475W), though this number is higher for galaxies 
near bright sources. We determine the uncertainty in the sky 
by dividing the width of the histogram by the square root of 
the number of boxes. We truncate the profiles at the radial 
distance where the uncertainty in the sky translates into a 
surface brightness uncertainty of 0.1 mag/arcsec 2 (the sky 
values subtracted by MultiDrizzle are usually of order 25 
counts in F814W and 30 counts in F475W, which is 2-3 
magnitudes fainter than for ground-based observations of a 
moonless sky). This means that in a typical case, without 
bright neighbours, the profiles are truncated at F814W (AB) 
surface brighness of 25.5 mag/arcsec 2 and 26.0 mag/arcsec 2 
in F475W. For a dwarf galaxy this means an extent of about 
2 to 3 effective radii, and a bit more for elliptical galaxies. 



4 MULTIPLE COMPONENT FITS WITH 
BAGATELLE 

For all of the galaxies in our sample we determined one- 
component structural parameters with a Bayesian code 
called bagatelle, which fits a parametrized light profile 
to a galaxy. For the light profile we assumed a Sersic profile, 
where we leave all parameters (Sersic index n, R c g, magni- 
tude, ellipticity, position angle, position) free, except the 
boxiness/disciness parameter, which is fixed to the value 
that generates perfect ellipses. A more elaborate description 
of the code can be found in a forthcoming paper (den Brok 
in prep.), but we summarize some of the characteristics of 
the code. 

Galaxies only rarely consist of one structural compo- 
nent. Even though a Sersic profile may be a good approxi- 
mation for the light profile of the outer galaxy, the centres 
of galaxies are known to have different profile shapes. For 
example, dwarf gala xies often show evidence for compact 
nucle ar star clusters (jGraham fc Guzman! 1200*31 ; ICote et al.l 
2006). Except for the fact that some profiles describe the 
surface brightness of a galaxy better, several authors have 
attached physical interpretation to profile shapes as well - 
for example, the separation of a galaxy in a Sersic profile 
and a star cluster, or a Sersic profile and a core (which may 
be scoured by black holes.) It is of course important that 
one uses a profile which describes the data best, however, in 
the case of 2d profile fitting, the reduced Xr is usually not a 
good measure of how well a profile fits. By adding more and 
more free parameters to a fit, the fit will become better, but 
probably not more physical. 

This problem can be solved by using Bayesian fitting, 
because with Bayesian fitting more complicated models (i.e., 
the ones with more parameters) are punished. We have 
therefore written a Bayesian code. The input for each model 
is a set of priors. The output a set of parameters together 
with a number called evidence, which tells you how much 
the prior volume has collapsed. The value of the evidence 
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Figure 2. 2-Component fit with BAGATELLE to a Coma dwarf 
galaxy. Shown are the different components (Sersic and point 
source). The Sersic-only fit is not shown, but is according to the 
difference in evidence (Aevidence = 151) extremely unlikely. 



can help us choose between models, because it ranks each 
model in order of likelihood. 

We have fitted our galaxies in Coma with two mod- 
els. One model consists of a single Sersic, the other of a 
single Sersic and a point source, since star clusters at the 
distance of Coma should not be resolved. The position of 
the point source is forced to coincide with the position of 
the rest of the galaxy, so that we do not end up fitting glob- 
ular clusters. In Figure [2] we show an example fit where the 
addition of a point source to the profile is a significant im- 
provement on the single Sersic fit. We have not fitted our 
galaxies with additional disc components. For the fainter 
galax ies in our sample, th is is in most cases not a problem 
(cf. iFerrarese et"al] ((2006 ) ) , however, we note that there are 
SO galaxies for which the Sersic index may be biased. We 
note that the magnitude of NGC4874 is underestimated be- 
cause of the small field of view of the ACS camera. As this 
galaxy is a cD galaxy and therefore has a disputable abso- 
lute magnitude anyway, we have decided not to correct for 
this. 



5 RESULTS 

5.1 Colour gradients 

In the following sections we present our results on colour gra- 
dients. We always use F814W (AB) and F475W (AB) mag- 
nitudes. Since B and Ic (in the Johnson- Cousins system) 



are more commonly used, we give here a few conversions: 

Ic - F814W(AB) = -0.38 

(from the WFPC2 Photometry Cookbook) 

B = F475W(AB) 

+ 0.329(F475W(AB) - F814W(AB)) + 0.097 

» F475W(AB) + 0.46 
B-Ia = 1.287(F475W(AB) -F814W(AB)) + 0.538 
(from lPrice et alj (120091 ) ) 

Ic = Ic(AB) — 0.342 
B = B(AB) + 0.163 

(from [Frei fc Gurml i|l994l )') In Figure [3] (top) we show the 
colour gradient as a function of host galaxy absolute F814W 
magnitude. Gradients are negative and steep in the large 
elliptical galaxies, the less luminous galaxies have shallower 
gradients which are often positive. The change from negative 
to positive gradients seems to happen just at the magnitude 
where there is also a well known change in the behaviour 
of structural parameters of galaxies, where galaxies fainter 
than My = — 18 (corresponding to w —19 in F814W) mag- 
nitudes follow a different size-lumin osity relation than more 
luminous galaxies (see for example iKormen 

M i|l985h . and 

for stellar populations IVader et alj ( 19881 )). At first glance 
this looks exactly what one would expect in a galaxy for- 
mation paradigm where galaxies form through monolithic 
collapse. 

The centres of galaxies are special environments. Dy- 
namical time scales are short here. Many galaxy centres 
show evidence for central massive objects such as supermas- 
sive black holes and nuclear star clusters. The formation of 
these objects is not well understood, and in particular their 
relation with the formation and evolution of the host galaxy 
is uncertain. If we want to understand what processes have 
shaped the stellar populations of the host galaxy to its cur- 
rent form, it is necessary to exclude the for this study ag- 
gravating contribution of any central AGN or star cluster. 

To do so, we carry our linear fits to the colour profile 
out again, but now excluding a the data within a radius 
of 3 pixels from the centre (75 pc at a distance of Coma). 
This is usually enough to exclude most of the light contri- 
bution of additional nuclear sources to the colour profile. 
However, inspection of the colour profile has shown that 3 
pixels is sometimes not sufficient ( see for example Fig[5J ■ For 
these galaxies, we manually altered the fitting range (max- 
imally excluding 10 pixels from the centre). A more com- 
plicated problem is determining gradients in large galaxies 
with structures such as bars, where an abrupt change in stel- 
lar populations between different morphological components 
can produce a strong colour gradient. These galaxies were 
also treated by hand, where we always used the outermost 
component to fit the gradient. 

In Figure [3] (bottom) we present the colour gradients 
for galaxies where we have excluded any central component 
as a function of host galaxy absolute magnitude. This re- 
analysis has had effect on elliptical galaxies as well as dwarf 
galaxies. Some elliptical galaxies have now less steep colour 
gradients, because we eliminated central reddening due to 
dust discs, whereas the to a large extent positive gradients 
that we observed for dwarf galaxies in Figure [3] (above) have 
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Figure 3. Top: Colour gradients of likely Coma members as a function of absolute magnitude. Bottom: 
Colour gradients of the same galaxies, now excluding the central parts from the fit. Red squares are 
compact galaxies from Price et al. (2009). The blue points required manual intervention. 



mostly disappeared and become slightly negative. There are 
two galaxies with strong positive gradients. One of these 
galaxies (CDMAil3005 . 684p275535 . 20) shows a central disk 
and spiral structure. This galaxy may be a spiral galaxy 
transforming to a dwarf galaxy, even though its projected 
distance is not too far away from the core of the cluster. 

Even without the positive gradients for dwarf galaxies, 
the trend is visible that fainter galaxies have smaller gradi- 
ents (a Spearman rank test confirms this with a correlation 
coefficient of 0.41 (excluding the compacts gives 0.45)). This 
is at least qualitatively consistent with a monolithic collapse 
scenario. In this scenario, after a first burst of star forma- 
tion, primordial gas is enriched with stellar ejecta until the 
onset of galactic winds prohibits any further star formation. 
Since the latter process is dependent on the local escape ve- 
locity, colour gradients are a function of host galaxy mass, 
for which we use absolute magnitude as a proxy in Figure [3] 

In a subsequent paper, we will investigate the relation- 
ship between the colour of the nuclear region with that of 
the main body of the galaxy. In the following part of this 
paper, we focus only on the gradient in the main body of 
the galaxy, i.e. excluding the central part. 



5.2 Differences in colour gradients as function of 
environment 

The environment in which a galaxy evolves can influence its 
colour gradient, since it can regulate the amount of star for- 
mation (by starving a galaxy from an external gas reservoir, 
by removing already present gas, by heating or ionizing cold 
gas or by compressing present gas due to the high-pressure 
in the ICM) and set the number of interactions with other 
galaxies, since harassment and mergers can alter the mor- 
phological appearance of a galaxy. One therefore expects to 
see differences between galaxies in the outskirts of the clus- 
ter and galaxies near the centre of the cluster. In Figure [4] 
we plot the colour gradients as a function of clustercentric 
distance. Although the results are not particularly sensitive 
to the way the distances are combined, we adopt as our 
indicator the harmonic mean of the distances to the two 
brightest galaxies, NGC4874 and NGC4889. It is difficult to 
find any environmental trend in Figure [4] Indeed the value 
of Spearman's rank correlation coefficient is 0.08 (0.06 when 
excluding the compacts) which means that there is no cor- 
relation at all. It is known from our previous analysis that 
more massive galaxies have stronger colour gradients. On the 
other hand, the clustering properties of massive galaxies are 
different from those of dwarf galaxies, in the sense that low- 
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Figure 4. Colour gradients of likely Coma members as function 
of harmonic mean distance towards NGC4874 and NGC4889. 
Red squares are the compact galaxies identified by Price et al. 
The size of the dot is related to the magnitude of the host galaxy 
as indicated in the figure. 



mass gas-poor galaxies are generally more strongly clustered 
than high-mass gas-poor galaxies. A better question to ask is 
perhaps if the colour gradients of galaxies with similar mass 
are somehow dependent on density. We therefore decided to 
divide out the gradient-magnitude relation for bright galax- 
ies, by fitting a line to the data of Figure [3] The residuals of 
this fit, i.e. the vertical offsets of the colour gradients from 
this line, are shown in Figure [5] as a function of distance. 
To check if there is any environmental dependence on the 
offsets, we have fitted a line through them, as a function of 
logarithmic projected distance. We find a slope consistent 
with zero, meaning that there is no evidence in our data 
that the offsets from the main relation are due to environ- 
ment. (Even when excluding the compact galaxies from the 
fit, the reduced % 2 remains high and the slope is consistent 
with zero). A Spearman correlation coefficient test shows no 
correlation (r = 0.06 or 0.04 without compacts.) 



6 COMPARISON WITH STRUCTURAL 
PROPERTIES 

The sequence of giant elliptical galaxies can be divided into 
two classes: the rapidly rotating low-luminosity ellipticals 
which are quite flattened, and their (generally) higher lu- 
minosity slow- or non-rotating counterparts. Despite their 
dynamical and morphological differences, elliptical galaxies 
are fitted exceptionally well by Sersics R~ surface brightness 
law, where n, known as the Sersic index correlates positively 
with galaxy magnitude and effective radius. 



6.1 Correlation with structural parameters 

In Figure[6]we show how the colour gradients correlate with 
the structural properties that we have derived with our 
2d Bayesian fitting code. We have used the results of the 
Sersic component from the best fitting model. If one only 
fits single-Sersic components to galaxies, one often finds a 



Figure 5. Offset of the colour gradients from a linear fit to 
the data in Fig. \3\ (bottom panel) plotted against the har- 
monic mean distance towards NGC4874 and NGC4889. Com- 
pact galaxies are red. The lines are fits through the sample (the 
red dashed line includes the compacts, the almost indistinguish- 
able black dashed-dotted line shows the fit after including the 
compacts.) 



too high Sersic index for dwarf galaxies. A lot of the be- 
haviour in this plot is expected. Smaller galaxies have less 
steep colour gradients (this was expected because fainter 
galaxies are generally smaller). Similar behaviour is observed 
in the panels where we plot the gradient against Sersic index 
and effective surface brightness. We have refitte d the galax- 
ies tha t were identified as compact ellipticals bv lPrice et al.l 
(2009). The offset in the compact sources is highest in the 
effective radius panel. Surprisingly, the one compact source 
with a positive gradient appears roughly consistent with the 
other sources. 



7 COMPARISON WITH LITERATURE DATA 

The colour of a stellar population is a function of both age 
and metallicity. Age, metallicity as well as other line index 
gradients have been studied. In this section we compare our 
colour gradients with Mgb gradients and metallicity gradi- 
ents, derived by various authors in the literature. The ad- 
vantage of comparing directly with index gradients is that 
we compare direct observables with each other. Compari- 
son with metallicity gradients involves a model to convert 
colour to metallicity, and requires assumptions on age gra- 
dients. For ages greater than 2 Gyr, and metallicities higher 
than [Fe/H] = —0.6, which is definitely the case for our com- 
parison galaxies, the conversion between colour and Mgb is 
completely linear, and independent of age and metallicity. 

However, for the interpretation of our colour gradients 
we also compare them to metallicity gradients. A compar- 
ison with previous work on colour gr adients is sh o wn in 
Fig. [71 The data s hown there are from IVader et ail (|l988T ) 
and lPeletier et al] (ll990T). The B- band magnitudes are taken 
from de Vaucouleurs et alj (|l99ll ) and converted to absolute 
magnitudes using the distances from the NASA/IPAC Ex- 
tragalactic Database (NED). There is clearly a big improve- 
ment in the range of magnitudes that is covered by our data. 



8 M. den Brok pf. al 



0.2 
0.1 
0.0 
-0.1 
•0.2 



0.2 
0.1 
0.0 
•0.1 
•0.2 



□ 


i 1 1 1 i 

• 




• 

i . . I 


1" 10" 

Effective radius 




1 , , ,1.1.1 
• 




^^^^ 

• 

i . . i .... 


1 4 10 

Sersic index (n) 


i i 

■ 


1 1 1 1 1 1 
• 


• 

• ■ 


• 

ii 



0.2 
0.1 
0.0 
•0.1 
•0.2 



18 19 20 21 22 23 24 25 26 27 

He (F814W) 



Figure 6. Colour gradients of likely Coma members, plotted against effective radius, Sersic index n and 
surface brightness at the effective radius. All quantities are calculated in the F814W band. The dashed- 
dotted line (V = —0.018 — 0.1372 log n) in the middle panel is a fit to the data, excluding gradients 
steeper than +0.15 and compact galaxies. 
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Figure 8. Colour gradients of likely coma members (black points) compared to metallicity and Mgb gradients from the literature (blue 
triangles). Upper panels: (left) Mgb data from Mehlert et al. (right) Mgb data from Sanchez-Blazquez et al. Lower panels: (left) Metallicity 
data from Sauron (Kuntschner et al. 2010) (right) Mgb data from Sauron (Kuntschner et al.) The Mehlert and Sanchez-Blazquez datasets 
plotted here consist entirely of Coma galaxies. However, the Sauron data has [Z/H] gradients compare better with our observations. 
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Figure 7. Comp arison of our colour grad ients in the outer part s 
with the work of lVader et~aH jl98g|) and lPeletier et al.l jl990l l. 



We do not confirm t he positive colour gradients for dwarf 
galaxies as found bv lVader et alj (Il988l ) , but our data are 
consistent with those of IPeletier et al.l l|l990t ). It is likely 
that previously correlations were not discovered because of 
this smaller magnitude range. However, it has been found 
by other authors that there is a turnover in p opulation gra- 
dient s around Mfsuw ~ —22. (for example, ISpolaor et al.l 
120091 . and references therein) . Unfortunately our our sample 
lacks galaxies in this magnitude regime. Since the galaxies 
with which we compare are quite large, seeing effects for 
the spectroscopic data should be negligible and exclusion of 
the central parts should in most cases not be necessary, as 
most of the bright galaxies do not have nuclear star clus- 
ters. For most of these data, only Johnson B magnitudes 
are available. We convert these to F814W magnitudes using 
B — F814W = 1.8 This is not entirely accurate, but the error 
which we make with this conversion should be less than 0.2 
mag for quiescent galaxies which are not affected by dust or 
recent star formation, and since we care mostly about global 
trends, we find this error acceptable. 
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7.1 Comparison with Mehlert et al. 2003 

iMehlert et all (|2003h fit line index gradients to data ob- 
tained by long-slit observations for a sample of 35 early-type 
galaxies in Coma. A benefit of comparing our data with 
this sample is, that we are sampling objects from a simi- 
lar environment, and any systematic error in the distance 
is the same for both samples . We used the galaxev soft- 
ware l|Bruzual &: Ch ariot 2003, henceforth BC03) to convert 
colours to an Mgb index. We measure the Mgb indices with 
the same index definition in galaxev, using the Padova 1994 
evolutionary tracks and a Salpeter IMF (lower and upper 
mass cut-offs at 0.1 and 100 Msun). From this analysis we 
find: 



AMgb(A) „ „ , , AF475W - F814W 

2_i — L = o 244 

Alog(r) Alog(r) 

Mehlert et al. define gradients in indices however as: 



A(Index) 



A log Index In (10) AIndex 



A log ft 



Index A log R 



(1) 



(2) 



Which means that there is an additional metallicity factor in 
the conversion from colour gradient to metallicity gradient, 
for which we use the signal-to-noise averaged central value 
from the same paper. A different way of doing this is fit- 
ting a relation between log(Mgb) and colour, but this gives 
essentially the same answer. Figure [8] shows a comparison 
between our derived gradients and the Mgb-derived colour 
gradients. The comparison shows that the scatter in these 
data is higher than in our data. Moreover, some gradients 
are positive, whereas the colour gradients in this magnitude 
regime are exclusively negative. 



7.2 Comparison with Sanchez-Blazquez et al. 
(2006) 

ISanchez-Blazauez et al.l l|2006t ) made an extensive study of 
metallicity gradients in early type galaxies, both for galax- 
ies in dense regions and field galaxies. Their sample includes 
also 21 members of the Coma cluster. Even though their full 
sample is much larger, again, to sustain homogeneity of the 
sample and avoid any additional uncertainties in absolute 
magnitudes due to distance uncertainties, we focus exclu- 
sively on their 21 Coma cluster galaxies. Because the step 
from index gradients to metallicity gradients requires a lot 
of assumptions on stellar models, we again just use the Mgb 
gradients measured by these authors. Because the gradient 
and index definition of Sanchez-Blazquez et al. is slightly dif- 
ferent from Mehlert et al., the conversion that we use here 
is just the simple: 



AMgb(mag) _ AF475W - F814W 

— U. J-00 - 



Alog(r) 



A log(r) 



(3) 



although in fact we use the inverse relation (i.e. we predict 
a colour from a given metallicity gradient). In the upper 
panel of Figure[S]we compare the predicted colour gradients 
from the Mgb line strengths to our observed colour gradi- 
ents. The converted gradients from Sanchez-Blazquez et al. 
are consistent with our observed gradients, except for two 
galaxies, which would have colour gradients steeper than 
—0.25 mag/dex. 



7.3 Comparison with SAURON 

iKuntschner et all |2010j) give a number of different line 
strength indices for early type galaxies in the SAURON 
sample, from which they determine metallicity and age gra- 
dients. Their data are obtained using integral field spec- 
troscopy, and have very high signal to noise. The SAURON 
sample consists of a magnitude limited sample of 48 early- 
type galaxies in the northern hemispere with t; ne j < 
3000km/s. 

Their fits are carried out in the same way as ours, by 
fitting a logarithmic gradient to the 2-dimensional metallic- 
ity maps, which were determined by fitting SSP models to 
different spatial bins. The SAURON sample contains galax- 
ies from nearby clusters and from the field, so a different 
density regime is probed by these observations. All galax- 
ies from the sample are well-studied nearby galaxies, so we 
assume that the distance uncertainty is small. 

Figure [8] (above) shows a comparison of the SAURON 
observations with our colour gradients. Amazingly, the data 
follows the same trends as ours: it shows a trend for gradients 
to become flatter for galaxies fainter than absolute F814W 
magnitude fainter than —20, the very shallow gradient for 
the very brightest galaxies (although we only have one data 
point) and also the scatter and magnitude of the gradients 
are completely consistent (if we exclude the one point with 
a very strong positive gradient.) 

We can also compare our colour gradients with the Mgb 
data from the same paper. Mgb is measured at two locations 
in each galaxy, at R e and R e /8. The fact that these gradi- 
ents are based on only two points, one of which is measured 
relatively close to the galaxy center, makes these gradients 
extremely sensitive to the presence of central stellar popu- 
lations which have formed much later than the galaxy halo, 
and hence the galaxies can show very strong gradients. In ad- 
dition to this, Mgb gradients can also become stronger if the 
a/Fe ratio increases toward the galaxy centers. The results 
in Figure [8] show that the Sauron data have a low intrin- 
sic scatter, which may be due to higher signal-to-noise than 
for example the Sanchez-Blazquez data. Also, the Sauron 
[Z/H] and Mgb gradients do not show the same trend. We 
conclude that at least in this magnitude regime our inter- 
pretation of colour gradients as metallicity gradients seems 
to be consistent with spectroscopic data. 



7.4 Comparison with Spolaor et al. 

ISpolaor et all (2009) use CMOS in long-slit mode to mea- 
sure spatially resolved line indices in 6 Fornax and 8 
Virgo cluster galaxies. They convert their Lick/IDS in- 
dices to meta ll icies and a ges using the SSP models from 
iThomas et all (|2003l . I 2004h . The metallicity gradients that 
they find drop from [Z/H] =— 0.3 to for fainter galaxies. The 
velocity dispersion-gradient relation found by these authors 
shows a surprisingly small scatter. We convert the metallic- 
ity gradient - B-magnitude data using BC03 to colour gra- 
dient - I-magnitude data. For this we use B(J) — I(AB)=1.8 
and AF475W - F814W = 0.27 A[Z/H] The results of this 
conversion are shown in Figure EH where we also overplo t 
the theoretical predictions from lKawata fc Gibson! (2003). 
Surprisingly, we find, if we compare the gradients with the 
results from Spolaor et al., that the results are consistent 
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Figure 9. Top: Colour gradients of likely Coma members and converted Spolaor et al. (2010) data. 
Bottom: Colour gradients after excluding the central part as function of absolute F814W magnitude. 
The blue points are the points from Spolaor et al., the dashed line is the theoretical prediction from 
Kawata & Gibson (2003) 



only if we use the full gradient, and not excluding the cen- 
tre. 



7.5 Comparison with Koleva et al. 

The comparison data in the previous subsections was mainly 
focussed on bright ellipticals. However, an important result 
of our work is that dwarf galaxies have slightly negative 
gradients. In this section we make a de tailed comparison be - 
tween our gradients and the data from I Koleva et al.l (|2009T ) , 
a paper studying metallicity gradients in dwarf galaxies. 
These authors obtained FORS1/2 long slit observations of a 
number of galaxies in the Fornax cluster, the Antlia cluster 
and the NGC5044 and NGC5898 groups. After rebinning to 
obtain sufficient S/N in the outer parts, they fit SSP models 
as a function of galactic radius to obtain ages and metallic- 
ities. 

Inspection of the data reveals that most age gradients 
are negligible, except for the centre - where dwarf and low 
mass elliptical galaxies usually host a star cluster - which 
is often, but not always, younger. We fit a gradient to the 
metallicity data, which we subsequently convert to a colour 
gradient using the BC03 models. Figure [TOl shows the gra- 
dients of the Fornax data compared to our own data. Re- 



markably, part of the galaxies seems to be consistent with 
our own gradients, whereas a group of points seems to be 
more consistent with the compact sources in Coma. 

For one Fornax galaxy, FCC136, ACS data is availabl e 
from the Fornax Cluster ACS Survey (| Jordan et all 120071 ) . 
in the F850LP and F475W passbands. The wavelength of 
the F850LP is a bit redder than our own ACS data, but this 
offset should be negligible for gradients. In Figure QT] we 
show a detailed comparison between the metallicity points 
and colour of FCC136. According to the data from Kol- 
eva, this galaxy has no age gradient, and put even more 
strongly, should even be a true SSP. Inspecting the colour 
profile shows that there is a nuclear star cluster present, 
which is much bluer than the rest of the galaxy, either be- 
cause of a lower metallicity, or because of a younger age. It 
is likely that the centre contains a much younger population 
that could not have been seen from the ground. 

Outside the centre, the metallicity profile does seem 
to follow the colour profile, with a small offset, which 
may be caused by uncertainties in the stellar popula- 
tion models. Foreground extinction in the direction of 
Fornax is small (a round E(B— V)=0.0155 according to 
ISchlegel et al.l dl998h) and hence E(g— z) should be around 
0.033 ( Sirianni et al.ll2005l ). This is a constant, but some- 
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from Koleva et al. converted to colour gradients. In black we 
show the data from the ACS Coma Cluster survey. 
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Figure 11. FORS metallicity data from Koleva et al. converted 
to colour data (red) with data from the ACS FORNAX survey 
in black. 



what uncertain colour offset, and we do not correct for it. 
(The small wiggle and subsequent flattening in the profile 
are due to use of a slightly too small PSF in the red band.) 
This plot shows also that one should be very careful in de- 
ciding which points to take into account for the fit. 



whole. Is there any ev idence for age gra d ients in 'normal' 
red-sequence galaxies? iKuntschner et al] (|2010h determine 
age gradients for galaxies in the Sauron sample and find that 
old (SSP age>8Gyr) galaxies are consistent with having no 
age gradients, but young galaxies often have a younger cen- 
tre. Sanchez-Blazquez et al. do find age gradients for their 
sample of galaxies in low density environments, but their 
results for galaxies in a high density environment (Coma) 
point at age gradients consistent with zero. This means that 
at least for the bright galaxies in our sample, the age gradi- 
ents are probably (on averag e) not important. However, we 
note that iRawle et all ()2010l ) warns that, even though they 
find on average zero age gradients, 40% of their sample is 
inconsistent with the absence of an age gradient. 
The situation may be very different in dwarf galaxies. 
IStinson"eT al. ( 2009) simulate age gradients in isolated dwarf 
spheroidal galaxies (dSphs) and find haloes which are sys- 
tematically older than the inner parts of the galaxy. Also 
oth er hydrodynamical simulations of dwarf ellipticals (see 
e.g. IValcke et al.l lioOc? ) show that star formation proceeds 
in bursts of gas that contracts towards the galaxy centre, 
with each subsequent starburst being more central, and even 
within one burst small age gradients are produced. The age 
gradient that is formed in this way is positive, and would 
lead therefore to a positive colour gradient. There has not 
been much observational work to test this with dwarf ellipti- 
cals. iKoleva et al. I (120091 1 measure gradients in bright dEs in 
the Fornax cluster, although they claim that most galaxies 
have positive age gradients, excluding the centre from the 
analysis the p r ofiles look remarkably flat. Also the results of 
ISpolaor et all (|2010h do not point at any age gradients. On 
the other hand, results for dSphs in the local g roup (see e.g. 
Battaglia etafl (|200fj )) and the M81 group (|Lianou et al.l 



2010T ) point at slightly older outer haloes of galaxies. 



If we assume that dwarf galaxies form outside-in (and 
this is where all the evidence from simulations points to), 
the fact that we see primarily negative colour gradients in 
dwarf galaxies tells us the following two things: dwarf galax- 
ies do have metallicity gradients and the contribution of the 
age gradient to the colour gradient is, although maybe not 
negligible, not as strong as the metallicity gradient. On the 
other hand, one of the two galaxies that has a strong posi- 
tive gradient may be a spiral galaxy transforming to a dwarf 
elliptical. The inverted gradient here may be due to recent 
star formation, so age gradients may be present if dwarfs are 
stripped spiral galaxies. 



8 DISCUSSION 

8.1 Metallicity gradients in dwarf galaxies 

One of the shortcomings of using broad-band data for study- 
ing gradients in galaxies is that one cannot separate age and 
metallicity contributions to the gradient. Both mechanisms 
are known to operate in galaxies. For example, M32 has flat 
colou r gradients in most colours (|Peletierll 19931 ; II auer et al.l 
1998) but shows both an age and metallicity gradient in 
long slit observations in the sense that the galaxy i s older 
and more metal poor on the outside l|Rose et al.| [2005). How- 
ever, M32 is classified as a compact elliptical and is probably 
not representitive for the class of red-sequence galaxies as a 



8.2 Comparison with theoretical models 

In the previous sections we have derived colour gradients for 
a homogeneous sample of Coma dwarf and elliptical galaxies. 
We have shown that the strength of the gradient is strongly 
dependent on the measurement method: including the cen- 
tre alters the gradients significantly. If we however exclude 
the centre, motivated by the fact that galaxy centres host 
cores, light excess, dust and star clusters, which, although 
being consequenses of evolution, are probably not generated 
by the initial collapse of the galaxy (assuming that a galaxy 
forms through a monolithic dissipational burst), the colour 
gradients stay slightly negative. The reduced gradients show 
a relation in which gradients become less steep towards faint 
magnitudes, but almost never become positive. This is at 



least qualitatively consistent with the old monolithic col- 
lapse scenario. Here, an elliptical galaxy forms as a large 
starburst from a primordial gas cloud. As not all the gas 
is used in the initial burst, the remaining gas sinks to the 
centre, whereas the stars stay where they are. On its way 
to the centre, the gas is enriched by stellar ejecta. In this 
way, a metallicity gradient is formed, with the stars near the 
galactic centre becoming more metal rich. If any age gradi- 
ent is present, it should be a positive one, but the metallicity 
gradient set up by the collapse is usually much stronger in 
the colour, an d there is no str ong observational evidence for 
age gradients. ICarlberd (|l984T) was one of the first to model 
dissipational collapse of an elliptical galaxy with an N-body 
code. In his simulation collisional gas cloud particles col- 
lide with each other to form collisionless star particles. The 
metallicity gradient of the simulated galaxies depends chiefly 
on the assumed ratio of the velocity of metal-enriched stel- 
lar ejecta and the velocity dispersion of the stars and has a 
value of VZ = —0.5 for typical elliptical galaxies. We note 
however that lowering this velocity dispersion yields gradi- 
ents that are more consistent with our observations, and we 
indeed confirm a trend that fainter galaxies have less nega- 
tive gradients. 

Since then models have improved and are b etter at re- 
producing observations. [ChIosi~&||Carrar3 (|2002T ) carried out 
N-body-tree-SPH simulations, thereby aiming at reproduc- 
ing the properties of dwarf galaxies. Their simulations start 
from a virialized proto-galaxy with a 1/r dark matter den- 
sity profile. One of their conclusions is that galactic winds 
are extremely important in blowing out metal rich gas from 
a galaxy and galactic winds of low mass galaxies are more 
efficient. From their Figure 6 we reconstruct metallicity gra- 
dients of approximately —0.8 for a 10 12 Mq galaxy and es- 
sentially zero for a 10 9 Mq galaxy. 

There have been many more attempts to model metal- 
licity grad ients, also in a cosmolog ical context. The simu- 
lations bv lKawata fc Gibson! (|2003T ) model elliptical galax- 
ies by following an isolated sphere with cold dark mat- 
ter density fluctuations superimposed. These simulations 
show the same trend from Chiosi & Carraro that less mas- 
sive galaxies have less steep gradients, which even go to 
zero for the least massive galaxies. Their simulation output 
for metallicity gradients does however contain only three 
different mass points for elliptical galaxies. Most recently, 
IPipino et al.l |2008, 2010) modelled metallicity gradients in 
elliptical galaxies. These authors reached the conclusion that 
the simple outside-in formation of ellipticals is insufficient to 
explain observables and requires more physics, such as de- 
tailed m odelling of gas flows . In st ark contrast with the mod- 
els from lKawata fc Gibson! (|2003), the Pipino et al. models 
show a variety of gradients, strongly dependent on the star 
formation efficiency. The galaxies with high SF efficiency 
have steeper gradients. Figure [12] compares our observed 
gradients with the theoretical models by Pipino et al. Their 
paper does not provide magnitudes of galaxies, so we make 
a simple con version to F814W ma gnitudes using the Milefl 
SSP models l|Vazdekis et al.ll2010T i where we assume a sin- 
gle age (10 Gyr) and metallicity ([Fe/H]=0) for all galaxies. 
Since we are only interested in the qualitative properties, we 

3 http://miles.iac.es 
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Eigenvalue 


Mag 


Eff Radius 


Sersic Index 


Colour 


3.17 


-0.54 


0.48 


0.52 


-0.46 


0.50 


0.14 


-0.64 


0.06 


-0.76 


0.25 


0.23 


0.54 


-0.67 


-0.46 


0.12 


0.80 


0.28 


0.53 


-0.05 



gradient 



Table 1. Results of a principal components analysis with struc- 
tural parameters. The first column gives the eigenvalue of each 
component, the other columns show the relative importance of 
the observables in this component. 

do not correct for the fact that smaller galaxies have on av- 
erage a lower metallicity. The range of metallicity gradients 
predicted by these models seems to coincide well with our 
observations. 



8.3 Structural parameters 

In section [5] we showed that colour gradients also correlated 
with structural quantities other than magnitude. It is impor- 
tant to know which of these quantities is the primary driver, 
since it might tell us how the build-up of the galaxy has 
proceeded. The strongest correlation that we found (even 
stronger than magnitude) was with Sersic index (Spearman 
rank coefficient 0.6)and as far as we are aware, the strong 
correlation of colour gradient with Sersic index is shown 
here for the first time. The Sersic in dex is a measure o f 
how centrally concentrated a galaxy is (|Truiillo et al.l[200lh . 
Although the exact physical processes that determine the 
Sersic index of a galaxy are somewhat clouded, we note that 
hig h Sersic indices can be produced by violent re laxation (see 
e.g. Ivan Albadal Il982l ; Irliorth fc Madsenl Il99ll ). As colour 
gradients are to first order a result of dissipational collapse, 
it is not so clear why they should correlate with a parame- 
ter resulting from a non-dissipational process. Other authors 
have also look e d at c orrelations with structural parameters. 
iTortora et all (|2010h claim to see a dip in colour gradient 
at Sersic index n=2. We do not confirm this with our data 
(the middle panel of Fig. [6]), but their data also contain 
late type galaxies, which may explain this difference. As a 
rule of thumb, quantities correlate better with velocity dis- 
persion than with any other structural quantity. All struc- 
tural quantities shown in this paper are know to correlate 
with one another at least in some way. It is therefore not 
so clear if any of these quantities is more fundamental. It 
is however possible that a combination of structural param- 
eters together with the colour gradients form a hyperplane 
with low scatter, which, when seen in projection against any 
individual structural parameter, has a consequently larger 
scatter. To test this, we performed a Principal Component 
Analysis of three of the structural parameters (magnitude, 
Sersic index n, effective radius) and colour gradient. The re- 
sults are shown in Table [1] The relative importance of each 
component is given by the eigenvalue of the component. The 
most dominant component is obviously the first component. 
The component divides its energy over the four variables in 
an almost equal way, indicating that there is no good evi- 
dence for a hyperplane, but more for a line. The second most 
important component shows that colour gradient is inversly 
related to galaxy size, but this component is not very strong. 
In Figure [6] we see that especially in the V-logr plot, the 
compact galaxies are outliers. The colours of these sources 
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Figure 12. Colour gradients of Coma galaxies (black dots and red squares) as a function of absolute F814W magnitude (and 
with an approximation of stellar mass) with the models from Pipino et al.(2010) superimposed. Simulated massive ellipticals, 
low-mass ellipticals and bulges are plotted with cyan, green and blue hexagons, respectively. The grey zone defines metallicity 
gradients between Pipino et al.'s maximum star formation efficiency-gradients and no gradient. 




jPrice et alj|2009h . colour gradients and surface brightness 
at the effective radius (Fig. [Gj are all similar to those of 
larger elliptical galaxies, and it appears that only their sizes 
are different. (The exception here is the smallest compact, 
which has apparently a positive gradient. However, given the 
extremely small size of this source, it is quite possible that 
this is the result of errors in the PSF.) Are the observations 
of our colour gradients consistent with a tidal threshing sce- 
nario? The statistics here are low, but it is tantalizing that 
the compacts seem to follow a similar V-log r relation, which 
is offset from the 'normal' relation by ~ 1 dex in r, but we 
do not see many galaxies in between the two sequences (as a 
matter of fact, there is an additional galaxy which seems to 
follow the compacts sequence on the large-end side). Obvi- 
ously, extrapolating these data back to find the progenitors 
of these galaxies leads to fairly large and bright galaxies, 
which would be difficult to strip tidally. It is therefore likely 
that, if these galaxies were tidally transformed, either ad- 
ditional central star bursts have happened, which made the 
galaxy metal rich in the centre, or the stripping process man- 
aged somehow to enhance an additional gradient. An addi- 
tional possibility is that we are looking at stripped bulges 
where still part of the disk is pre sent (note that there is ev i- 
dence for disc signatures in M32 (|Grahamll200l ; ISethll2010h ) 
In Figure [T3] we show how the colour gradient-magnitude 
relation depends on Sersic index, by giving the points a 
different colour for different Sersic indices. This plot gives 



the impression that in a given magnitude bin, galaxies with 
higher Sersic index have a steeper colour gradient. We warn 
however that we have not carried out a bulge-disc decom- 
position for our galaxies and therefore the interpretation of 
the Sersic index may be different for spheroids with disc 
than for spheroids without disc. One of the problems of our 
sample is that it is a mixed bag of objects. The presence 
of multiple morphological components in a galaxy impedes 
a clean interpretation in terms of formation scenarios for 
elliptical galaxies. As a check, we plot in Fig. [14] a set of 
galaxies for w hich morphologies have been determined by 
iDresslerl JT980) . Unfortunately, these morphologies are only 
available for the brightest galaxies. In this figure, we see that 
in general elliptical galaxies have steeper gradients than SO 
galaxies. We can speculate that rotational support has pre- 
vented gas from flowing to the centre and therefore thwarted 
the formation of a strong metallicity gradient. The outlier 
in this plot is NGC4876, a galaxy with overal quite a flat 
gradient but a very red centre due to a dust ring. 

8.4 The influence of the cluster on gradients 

If the colour gradient of a galaxy is influenced by the star 
formation efficiency (which in turn is dependent on tidal 
interactions, pressure in the ICM etc.) we should expect 
some dependency on clustercentric radius (or at least for 
the most massive galaxies.) However, since the virial radius 
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Figure 13. Colour gradients of Coma galaxies of absolute mag- 
nitude. The colouring denotes the Sersic index of the host 
galaxy. Interestingly, there is the suggestion that in a given 
magnitude bin, galaxies with higher Sersic index have a steeper 
colour gradient. Compact galaxies are marked with a red cross. 
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(2009) shows that some of the galaxies in the Fornax cluster 
(which in terms of density resembles more a rich group than 
a cluster) have much stronger metallicity gradients than our 
Coma dwarf galaxies a t the same magnitude. T his is consis- 
tent with the result of lLa Barbera et al.l l|2005l ) who found 
that the steepness of colour gradients is dependent on the 
richness of the cluster, although the galaxies observed in 
that paper are much brighter. It is possible that there are 
different formation scenarios for dEs, and that therefore the 
distribution of colour gradients between our Coma sample 
and Koleva et al.'s Fornax sample is slightly different. 

8.5 On the necessity of high resolution data 

Our new result which shows that dwarf galaxies have pri- 
marily flat or negative gradients is in contr ast with what 
other authors have found in the past (see e.g. Ivan Zee et al.l 
l2004h and we have shown that this is primarily due to the 
presence of nuclei and dust discs in the centres of galax- 
ies. High resolution observations have shown that almost all 
dwarf galaxies have central stellar nuclei, which have dif- 
ferent colours than the rest of the galaxy. The fact that 
Ivan Zee et al.l l|2004h detect positive colour gradients in non- 
nucleated dwarf galaxies has probably more to do with a blue 
nucleus going unnoticed than that the main component of 
the galaxy really has a positive colour gradient. 

It is well known that galaxies can contain more than 
one component. For elliptical galaxies and dEs, we can often 
fit a light profile with a single Sersic profile, if we exclude 
the centre, so that we are not plagued by the presence of 
cores, light excess or nuclear star clusters. However, our data 
indicates that the colour profile is not always a straight line, 
but can show breaks. These coincide with the presence of 
other morphological components, for example, the presence 
of a nuclear star cluster or a central nuclear disc. Separating 
the components, and individually studying the photometric 
properties of them may lead to a better understanding of 
the formation of these galaxies. 



Figure 14. Colour gradients of Coma galaxies for which mor- 
phologies are available (from Dressier (1980)), as a function 
of absolute magnitude. The colouring denotes the Sersic index 
of the host galaxy. Ellipticals have steeper gradients than SO 
galaxies. The outlier among the ellipticals is NGC4876. 

of the cluster is 2.9h^Q Mpc ()Lokas fc M amon 20031). we are 
only probing the very densest part of the cluster, so a sim- 
ple explanation may be that the baseline in densities over 
which we are comparing our gradients is simply not large 
enough. In addition to this, the fields around NGC 4839 are 
comparable to the density of a rich group, and it may be 
the case that the region is sufficiently dense such that star 
formation has been switched of here just after a first burst. 
Besides this, we only have a small number of galaxies which 
we can study. To really investigate the effect of environment, 
a comparison with galaxies in a much lower density would 
be the next step. Theoretically one would at least expect 
that tidal processes which change the orbits of stars to more 
radial orbits, at least slightly flatten the colou r gradient. A 
comparison with the metallicity gradients from I Koleva et all 



9 SUMMARY AND CONCLUSIONS 

We have derived colour gradients for a sample of confirmed 
or very likely Coma cluster members. We show that the gra- 
dient which one measures is strongly dependent on whether 
or not one includes the centre of a galaxy. Large ellipticals 
often have dusty discs in their centre, whereas lower mass el- 
liptical or dwarf galaxies have nuclear star cluster, of which 
the colour is different from the rest of the galaxy. This shows 
that for accurate measurement, one really requires high res- 
olution observations. 

Our findings indicate that colour gradients in general 
correlate with galaxy luminosity or mass, in the sense that 
high-mass galaxies have stronger gradients than dwarfs. 
However, in general, colour gradients in the main body of 
the galaxy do not become positive, as has been claimed in 
the literature. The colour gradients seem to correlate equally 
well with other structural parameters, such as effective ra- 
dius, effective surface brightness and in particular Sersic in- 
dex. Compact galaxies clearly stand out by having steeper 
gradients than can be expected on the basis of their mag- 
nitudes. For a subset of galaxies of which the morphologies 
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are known, we find that SO galaxies have less steep gradients 
than elliptical galaxies. 

We have compared our colour gradients with Mgb gra- 
dients and metallicity gradients from the literature, by con- 
verting these gradients to colour gradients. The comparison 
between colour gradients and Mgb gradients is consistent, 
but the Mgb gradients show generally much larger error- 
bars. Recently, Spolaor et al. measured an extremely tight 
relation between metallicity gradient and mass. A compari- 
son of these data with our data shows results consistent with 
our colour gradients if we include the centre in our fit, but 
not if we measure the gradient in the outer parts. Since time 
scales are much longer in the outer parts, and later bursts of 
star formation primarily influence the centre, the gradient 
in the outer parts should be a better probe for metallicity 
gradients from a monolithic collapse scenario. The metal- 
licity gradients from the SAURON data are however fully 
consistent with our data. 

Simulations still have tr ouble reprodu c ing th e observa- 
tions. Although the work bv lPipino et al.l Q2010) is able to 
reproduce the scatter, it apparently misses the trend ob- 
served with magnitude. If dwarf galaxies have strong posi- 
tive age gradients, the results are consistent. 

Given the poor reproduction of gradients by simula- 
tions, we can only speculate what our results mean in 
terms of galaxy formation models. The observed trends 
with galaxy magnitude imply that somehow the mass of 
the galaxy, or the potential well are important in shaping 
the gradient. However, given the scatter in gradients at a 
given magnitude, this cannot be the only important process. 
Another clue comes from the Sersic index-gradient relation. 
Galaxies with nearly exponential profiles have flatter gradi- 
ents than galaxies with higher Sersic indices. Higher Sersic 
indices are thought to be the results of processes involving 
violent relaxation. Since this is a non-dissipational process 
by itself, the strong gradients may point at a history of moist 
or wet mergers. 
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195, 


17021 


27, 


,99662 


13, 


.97 





0, 


,0236 


-0.045±0.004 





,31 


0.15 


COMAil3041 


.192p28242.38 


195, 


17162 


28, 


04510 


16, 


.76 








,0287 


-0.047±0.003 


0. 


,22 


0.15 


COMAil3042 


.519p28325.45 


195, 


17719 


28, 


,05707 


18. 


.75 








,0195 


-0.153±0.013 


0. 


.05 


3.36 


COMAil3042 


.833p275746.98 


195, 


17847 


27, 


,96303 


14, 


.11 








.0280 


-0.060±0.001 





,20 


2.11 


COMAH3042 


.889p28313.73 


195, 


17868 


28, 


05382 


18. 


.02 








.0226 


-0.051±0.010 





33 


0.71 


COMAil3044 


.127p28215.37 


195, 


18385 


28, 


,03760 


18 


.46 








.0292 


-0.005±0.020 





,17 


2.01 


COMAil3044 


.634p28602.30 


195, 


.18599 


28, 


10064 


14 


,89 








.0219 


-0.089±0.000 





,48 


0.14 


COMAil3048 


,045p28557.42 


195, 


20018 


28, 


,09928 


18 


,08 








.0208 


-0.012±0.051 


0. 


.10 


0.98 


COMAil3048 


.646p28526.74 


195, 


20268 


28, 


,09075 


13. 


,43 





0, 


.0230 


-0.117±0.007 





,25 


1.71 



Table Al: For each galaxy, Coma ID (see Paper II), right ascension 
and declination are given, together with Sersic magnitude, membership 
ratio, redshift, colour gradient in the outer parts, ellipticity used in the 
Galphot fit and the reduced \ 2 from the fit of the gradient to the outer 
part of the colour profile. 
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COMAil3017.020p28350.10 

COMAil3018.349p28333.28 

COMAil3011.807p28502.62 

COMAil3039.113p28035.46 

COMAil3042.508p28324.94 

COMAil3037.589p28059.94 

COMAil25945.531p28312.61 

COMAil3013.398p28311.81 

COMAil25914.447p28216.87 

COMAil3042.753p275816.88 

COMAil3017.685p275718.92 

COMAil3018.094p275723.51 

COMAil3024.849p275921.86 

COMAil25959.071p275841.26 

COMAil25944.776p275807.13 

COMAil25930.832p275810.32 

COMAil25942.931p275954.39 

COMAil25956.707p275548.62 

COMAil25950.703p275514.94 

COMAil25944.825p275536.87 

COMAil25923.730p275042.66 

COMAil25814.969p272744.81 

COMAil25833.136p272151.77 

COMAil25831.666p272342.04 

COMAil25837.251p271035.59 

COMAil25825.308p271200.04 

COMAil25634.648p271339.03 

COMAil25623.788p271402.30 

COMAil25638.099p271304.09 

COMAil25857.453p274706.28 

COMAil25854.632p274742.16 

COMAil25842.637p274537.86 

COMAil25844.579p274458.30 

COMAil25845.297p274650.75 



late type, bar flattens gradient 

gradient in sky, maybe positive gradient 

too low S/N 

irregular galaxy 

in edge of a frame 

too faint (F814W = 19.9) 

in edge of a frame 

completely edge on 

too faint (F814W = 19.8) 

companions + bar + dust 

N4898a (overlapping) 

N4898b 

in low S/N gap between chips 

too low S/N 

too low S/N 

strong gradient in sky 

sky very uncertain 

overlapping, spiral arms 

strong gradient in sky 

too faint (F814W=19.99) 

in edge 

late type 

late type 

late type 

late type 

late type 

late type 

late type 

late type 

in edge 

in edge 

late type 

overlapping 

overlapping 



Table A2. Galaxies excluded from our sample. 
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Galaxy 


Aage 


A[M/H] 


Ills 1 


AMgb 


Acol 


NGC4673 


-0.1593 ±0.096 


-0.1049 ±0.1239 


13 


7 


-0.023 ±0.003 


-0.15 


±0 


019 


NGC4692 


0.163 ±0.1103 


-0.3817 ±0.1285 


14 





-0.014 ±0.004 


-0.09 


±0 


026 


NGC4839 


0.28 ±0.1708 


-0.3535 ±0.154 


13 


(2) 


-0.016 ±0.02 


-0.10 


±0 


129 


NGC4842A 


-0.2147 ±0.2617 


0.0637 ±0.2364 


14 


9 










_ 


NGC4864 


0.49 ±0.1762 


-0.5792 ±0.1618 


14 


8 


-0.016 ±0.009 


-0.10 


±0 


058 


NGC4865 


0.0672 ±0.1232 


-0.2587 ±0.1416 


14 


6 


-0.026 ±0.01 


-0.17 


±0 


065 


NGC4874 


0.5845 ±0.1029 


-0.7181 ±0.1382 


13 


7 


-0.044 ±0.01 


-0.28 


±0 


065 


NGC4875 


0.1777 ±0.139 


-0.3231 ±0.2141 


15 


6 (3) 










_ 


NGC4889 


0.2578 ±0.0711 


-0.4363 ±0.112 


13 





-0.012 


±0.005 


-0.077 


±0 


032 


NGC4908 


0.0131 ±0.2231 


-0.2354 ±0.2812 


14 


9 


-0.004 


±0.005 


-0.026 


±0 


032 


IC832 


-0.1432 ±0.1404 


0.0902 ±0.1483 


15 


(3) 


-0.003 


±0.004 


-0.019 


±0 


026 


IC3957 


0.2298 ±0.2385 


-0.4788 ±0.2636 


15 


6 (3) 


-0.006 


±0.009 


-0.039 


±0 


058 


IC3959 


0.0731 ±0.1344 


-0.1642 ±0.1655 


15 


2 


-0.009 


±0.004 


-0.058 


±0 


026 


IC3963 


0.336 ±0.134 


-0.5703 ±0.1788 


15 


7(3) 


-0.011 


±0.008 


-0.071 


±0 


052 


IC3973 


0.4122 ±0.1776 


-0.3849 ±0.1643 


15 


2 


-0.005 


±0.004 


-0.0323 


±0 


026 


IC4042 


-0.1029 ±0.4661 


-0.3628 ±0.6512 


15 


4 


-0.05 


±0.007 


-0.323 


±0 


045 


IC4051 


-0.0582 ±0.0852 


-0.3693 ±0.1261 


14 


8 


-0.03 


±0.014 


-0.194 


±0 


090 


CGCG159-41 


0.1629 ±0.5742 


-0.0679 ±0.6968 


15 


5 


-0.022 


±0.005 


-0.142 


±0 


032 


CGCG159-43 


0.5593 ±0.1614 


-0.6726 ±0.1936 


15 


3 


-0.014 


±0.004 


-0.090 


±0 


026 


CGCG159-83 


-0.4099 ±0.131 


0.0939 ±0.1566 


14 


9 


-0.023 


±0.008 


-0.148 


±0 


052 


CGCG159-89 


-0.0815 ±0.141 


-0.3979 ±0.1969 


14 


8 


-0.016 


±0.004 


-0.103 


±0 


026 



Table A3. Data from Sanchez-Blazquez et al. extended with own 
data. Columns show age, metallicity and Mg b gradien ts, w i th for - 
mal errors. B-band magnitudes are all from f DlFalco et al.l lll999T) 
except (2) I Gil de Paz et al] d2007h and (3) Ide Vaucouleurs et all 
dl99lh 
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Figure Al. Colour profiles for a subset of our galaxies, sorted by magnitude in anticlockwise direction. Brightest galaxy 
here is CDMAil3022 . 156p28249 . 08 (F814W(AB)=13.77) and faintest galaxy is C0MAH25902 . 433p28021 . 36 (F814W(AB)=18.93). 
COMAi 13022 . 156p28249 . 08 is a barred galaxy (which is why the errors increase and decrease over the profile) and has a cen- 
tral dust disk.Also C0MAi25943 . 724p275940 . 89 has a central dust disk, which is seen face on. COMAil25902 .433p28021 . 36 and 
COMAil25635.495p271490.37 are typical examples of dwarf galaxies, whereas C0MAH3007 . 123p275551 . 49 is an example of a faint system 
with a steep gradient and high Sersic index. The profiles of the full sample can be found in Appendix B of the electronic edition of this 
paper. 



